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ABSTRACT 

We present VLBA observations and a statistical analysis of 5 GHz VLBI 
polarimetry data from 177 sources in the Caltech-Jodrell Bank flat-spectrum 
(CJF) survey. The CJF survey, a complete, fiux-density-limited sample of 293 
extragalactic radio sources, gives us the unique opportunity to compare a broad 
range of source properties for quasars, galaxies and BL Lacertae objects. We 
focus primarily on jet properties, specifically the correlation between the jet axis 
angle and the polarization angle in the core and jet. A strong correlation is 
found for the electric vector polarization angle in the cores of quasars to be 
perpendicular to the jet axis. Contrary to previous claims, no correlation is 
found between the jet polarization angle and the jet axis in either quasars or BL 
Lac objects. With this large, homogeneous sample we are also able to investigate 
cosmological issues and AGN evolution. 

Subject headings: galaxies: active - galaxies: jets - galaxies: nuclei - radio 
continuum: galaxies; polarimetry 



1. Introduction 

Studies of the polarization properties of jets on the milliarcsecond scale can yield vital 
insights about their formation, collimation, and propagation (e.g. Gomez et al. 2000; Meier 
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2001). Based on results from 25 BL Lacs, Gabuzda, Pushkarev & Cawthorne (2000) claim a 
tendency for the electric vector polarization angle to lie along the jet axis for these objects, 
with about 30 percent showing a perpendicular orientation. This has been interpreted in 
terms of oblique and transverse shocks (Gabuzda, Pushkarev & Cawthorne 2000) or as the 
result of a helical magnetic field confining the jet (Gabuzda 2002). 

Although many individual sources and even large surveys have been undertaken with 
VLBl in total intensity, the number of sources studied with VLBl polarimetry is far less. In 
the early years of VLBl this was due in part to hardware restrictions (e.g., feeds with only 
a single available polarization), but also due to software restrictions that made calibration 
difficult. The VLBl group at Brandeis pioneered efforts to develop techniques and software 
to make VLBl polarimetry more tractable (Roberts et al. 1984; Roberts, Brown, & Wardle 
1991), and this culminated in a study of the mas polarization structure of 24 sources from 
the Pearson- Readhead sample (Cawthorne et al. 1993). But it wasn't until the advent of the 
VLBA, with homogeneous, low-leakage antennas and feeds, that VLB! polarimetry became 
widely accessible. Full polarization calibration techniques were implemented in AlPS in the 
1990s and in 1999 NRAO began providing regular monitoring^ of a suite of calibrators that 
can be used for the absolute polarization angle calibration (Taylor & Myers 2000). 

The Caltech-Jodrell Bank fiat-spectrum (CJF) survey is a large complete, fiux-limited 
sample of sources imaged with VLBl, consisting of 293 sources (Taylor et al. 1996). This 
survey has been used to place constraints on gravitational lensing by 10^ Mq black holes 
(Wilkinson ct al. 2001), to study the unbcamcd synchrotron luminosities of relativistic jets 
(Lister & Marscher 1997), and recently to study the relation between the linear size of the 
VLBl jet and the black hole mass (Cao & Jiang 2002). 

As part of a proper-motion study of compact objects, all the sources in the CJF sample 
were observed at least 3 times at intervals of approximately 2 years (Britzen et al. 1998). 
The early epochs used Global VLBl observations, but starting in 1995 the VLBA was used 
exclusively since it could provide better {u, v) coverage and consequently improved image 
fidehty. The first observing sessions with the VLBA were hmited in bandwidth (8 MHz) and 
used only a single polarization, but as the VLBA matured and its capabilities increased, we 
began in 1998 to observe with 16 MHz bandwidth in full polarization mode. The last three 
epochs of CJF were observed in this mode and provide full VLBl polarimetry of 177 sources. 

Here we present a statistical analysis of the properties of those 177 sources in CJF 
observed with full polarimetry. The sample is described in detail in §2, and the calibration 
procedures are described in §3. We present results in §4, and discuss comparisons of source 
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properties, and evolution of those properties with redshift in §5. 
We assume Hq = 75 km s~^ Mpc~^ and go=0.5 throughout. 

2. Sample Selection 

We present 5 GHz Very Long Basehne Interferometric (VLBI) polarimetry observations 
of 177 sources in the Caltcch-Jodrcll Bank flat-spectrum (CJF) survey, a complete flux- 
density-limited sample of 293 flat-spectrum radio sources. The CJF sample requires that 
the flux density at 4850 MHz be at least 350 mJy, (54850 > 350 mJy), and that the spectral 
index at 4850 MHz and 1400 MHz be at least -0.5 {afH^ > -0.5). Additionally, the sample 
is hmited by B1950 dechnations and galactic latitudes such that S > 35° and |&| > 10°. A 
list of the 293 extragalactic radio sources and a more complete description of the selection 
criteria for the CJF sample can be found in Taylor et al. (1996). The CJF sample has 
been completely observed with VLBI as part of the PR (Pearson & Readhead 1988), CJl 
(Polatidis et al. 1995) and CJ2 (Taylor et al. 1994) surveys. 

The 177 sources presented here were observed as part of a proper motion study (Britzen 
et al. 1998) so that polarimetric data in multiple epochs are available for 66 of the 177 
sources. Where data from multiple epochs was available, we preferred the 2000.958 obser- 
vations over the 1998.122 observations, and used the 1999.890 observations, with somewhat 
poorer u,v coverage, only when other epochs were unavailable. All plots shown here include 
just one epoch of each of the 177 sources. In the future it may be interesting to study 
variability in polarized intensity and polarization angle using those 66 sources with multiple 
epochs available, however no effort has been made toward this goal thus far. 

Using optical classiflcations from the literature (Henstock et al. 1995; Vermeulen, Taylor, 
Readhead, & Browne 1996), from the NASA Extragalactic Database, and from unpublished 
observations by the CJ collaboration, we find that the final sample of 177 sources consists of 
30 galaxies, 106 quasars, 20 BL Lacertae objects and 21 others, including empty fields, red 
and blue objects and other optically unresolved objects. We note here a few words of caution 
regarding the above classiflcations in that low-luminosity BL Lacs can be missidentified as 
quasars or radio galaxies (Marcha & Browne 1995), and variability in the emission lines 
of BL Lacs and quasars can cause some movement between these classes (Vermeulen et al. 
1995). Redshifts have been measured for 147 of the sources presented here. 
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3. Observations 
3.1. Calibration 

The observations were carried out on 2000 December 16 and 17 (2000.958), 1999 Novem- 
ber 21, 23 and 26 (1999.890), and 1998 February 8, 12, 13, 20 and 21 (1998.122). These 
three epochs provided about 34, 72 and 120 hours of data, respectively, for a total of 226 
hours of observations. The 2000.958 and 1998.122 observations used the 10 element Very 
Long Baseline Array"^ (VLBA) while the 1999.890 observations used only 8 of the VLBA 
antennas. The Saint Croix, Virgin Islands antenna was lost due to hurricane Lenny, and 
the North Liberty, Iowa antenna was lost due to a power hardware problem. Because we 
did not use phase referencing, we failed to image 6 sources with total intensity peaks of less 
than 55 mJy. Of these, four sources were identified with galaxies, and two sources had no 
optical identification. Another source, 0954+556, was omitted due to a position error of 30 
arcseconds. Finally, the gravitationally lensed object 0218+357 was also omitted. These 8 
sources that we did not image have not been included in any of the distributions or analysis, 
and are not included in the 177 source count. 

We observed at 4995 MHz with a total bandwidth of 16 MHz. Right- and left-circular 
polarizations were recorded using 1 bit samphng at 64 Mbits/sec. Two intermediate frequen- 
cies (IFs) of 8 MHz each were used for right- and left- polarization. Amplitude calibration 
for each antenna was derived from measurements of the antenna gain and system tempera- 
tures during each run. Global fringe fitting was performed using the AlPS task FRING, an 
implementation of the Schwab & Cotton (1983) algorithm. The fringe fitting was performed 
on each IF and polarization independently using a solution interval of 2 minutes, and a point 
source model was assumed. Next, a short segment of the cross hand data from the strongly 
polarized cahbrator 3C 279 was fringe fitted in order to determine the right-left delay differ- 
ence, and the correction obtained was applied to the rest of the data. Once delay and rate 
solutions were applied the data were averaged in frequency over 8 MHz. All data were then 
averaged together over a 30 second time interval. Imaging, editing and self-calibration were 
performed in DiFMAP (Shepherd, Pearson & Taylor 1995) and AIPS. 

Determination of the leakage terms was performed with the AIPS task LPCAL using 
a strong cahbrator with simple source structure such as 0923+392 or 0716+714. Given the 
nature of the survey, strong sources with good parallactic angle coverage were abundant 
in every observing session. The absolute electric vector polarization angle (EVPA) calibra- 
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tion was performed using various well-known target sources such as 0923+392, 1803+784, 
and 2200+420 (BL Lac). These observations were compared to contemporaneous short 
VLA observations and VLBI polarimetry reported in the literature, or available from the 
VLA/VLBA Polarization Calibration Page (Taylor & Myers 2000). We also used compo- 
nent C4 of 3C279, which was fairly stable during this time period (Taylor 2000; Zavala & 
Taylor 2001). The absolute uncertainty in the EVPA calibration is ~4°, with the dominant 
source of error being the variabihty of the EVPAs of the calibrator sources (Taylor & Myers 
2000). 

3.2. Data Analysis 

We automated the imaging process by writing a script that ran the AlPS task IMAGR 
on the self-calibrated data for every source, creating Stokes I, Q and U maps. Using another 
script, we then ran the AIPS task COMB on each of the sources to produce polarized intensity 
maps, polarized intensity noise maps, polarization angle maps and polarization angle noise 
maps. Some sample images covering a range of morphologies are shown in Fig. 1. These 
FITS files were brought into IDL where we used programs to objectively and uniformly 
determine the properties of each source. 

3.3. Source Property Definitions 

The RMS noise levels were calculated for the total intensity and linear polarizations 
by taking the standard deviations of four corner sections of the image and averaging the 
two middle values. We calculated the integrated flux in Stokes I, Q and U by summing 
the flux in a particular region that completely enclosed each source. The algorithm used to 
systematically define the region boundaries creates a rectangular box around those pixels 
with values at least 15 times the RMS noise level. Due to image artifacts this high cut- 
off proved necessary to correctly enclose many sources. A border of width three times the 
FWHM of the convolving beam was added to the initial rectangular area to define the final 
region to be used in the summing. As a check we compared the integrated fiux calculated 
in IDL to the sum of the clean components found in the AIPS task IMAGR. The total 
integrated fiux measurements of all 177 sources were found to agree with that calculated by 
the AIPS clean algorithm to within 8%. 

The core of each source was assumed to be located at the brightest peak in the total 
intensity map when no other data were available for study. However, in the case of compact 
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symmetric objects (CSOs) or sources with complex jet structures we took advantage of 15 
GHz images (Taylor, Readhead & Pearson 1996) (and unpublished 15 GHz observations) 
and published motion studies (Taylor & Vermeulen 1997) to more accurately determine core 
positions. Jet component positions were found by subtracting from the initial total intensity 
image a Gaussian model centered at the peak flux position and scaled so that the maximum 
value of the Gaussian equaled the peak value in the map. As the appropriate Gaussian 
is subtracted from the previous image, the locations of each new peak define the locations 
of consecutive jet components. We add the restriction that the total intensity at the jet 
component position found must be greater than 9 times the RMS noise level, where the 
noise is calculated as discussed above. At most 8 jet component positions were required to 
describe a source. The position angle of the elliptical Gaussian model used in the subtraction 
was matched to the position angle of the convolving beam, and both axes of the Gaussian 
were given a width of 1.5(T, where a = ^^JiM. and describes the width of the convolving 
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beam. Due to slight differences in the editing process that changed the u,v coverage for 
different Stokes parameters, the beam sizes and position angles fit by the AlPS task IMAGR 
differed slightly in I, Q and U. For consistency we forced all beams to have the dimensions 
and orientation of the Stokes U beam by setting the beam parameters BMAJ, BMIN and 
BPA in AIPS for the Stokes Q and I images. 

The jet properties and source morphologies based on total intensities were easily char- 
acterized after each jet component position was defined. We found 34 sources with no jet 
components which we call naked cores and 134 sources with at least one jet component which 
we call core-jet sources. In addition, we found 9 sources which we classified as Compact Sym- 
metric Objects (CSOs) or CSO candidates based on their symmetric structures. CSOs are 
a recently identified class of radio sources smaller than 1 kpc in size with emission on both 
sides of the central engine, and are thought to be very young objects (~1000 y, Readhead 
et al. (1996); Owsianik & Conway (1988)). Eight of these sources were previously identified 
as CSOs or CSO candidates (Peck & Taylor 2000), and the remaining source (0402+379) is 
currently under study. 

We recorded the jet axis angle to be the angle defined by the core and the closest jet 
component, where North is 0°. The only exception to this is the CSO candidate 0402+379 
which has a component close to the core but seemingly unrelated to the jet. The jet axis 
angle for this source was measured by eye. The jet length is defined in most cases as the 
distance from the core to the farthest jet component from the core, irrespective of jet bend. 
It is quoted in milliarcseconds and is uncorrected for angular distances. In §5 we model our 
angular size calculation after the description given in Kellermann (1993), creating a contour 
line at 2% of the peak in the total intensity map. The maximum distance from the contour 
line to the core position defines the angular size. To avoid measuring distances to spurious 
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peaks above the 2% level, we confined the contour to be within the same rectangular area 
that was previously determined to completely enclose the source. 

We define integrated core fiux densities as the pixel value at the core component, where 
the pixel values have units of Jy beam~^ and we assume the core itself to be unresolved. 
Due to complex jet structure, rather than sum the pixel values for the jet components we 
define the integrated jet fiux as the difference between the integrated total source fiux and 
the core fiux. 

Due to leakage between orthogonal polarization feeds in the antennas, we define detected, 
polarized cores or jet components as only those components in which the polarized intensity 
is greater than 0.2% of the peak in the total intensity image {p > 0.002 Ipcak)- We require 
detected polarizations to be above the Sa level, where cr is described by the RMS noise in 
the Stokes Q and U maps by the relation a = (0.21(7^ + 0.21 (Jq)^/^ according to Rayleigh 
statistics. In the case that the core has undetected polarization, we quote mcore as an upper 
limit equal to the larger of the two detection restrictions (mcore < max(0.002 Ipeak, Scr)). 
When studying correlations between jet axis angle and polarization angle, we restrict our 
analysis to those components with cr^ < 20°, where is determined by the noise image 
created by the MPS task COMB. 

4. Results 

Measured properties for each source are presented in Table 1. Additional information, 
including measured jet properties, is available on the world wide web^. Images are available 
upon request to the authors. Below we present a statistical analysis of source properties. 

4.1. Core Fraction 

From Fig. 2 it is clear that the core fractions of galaxies come from a different parent 
distribution than that of quasars or BL Lacertae objects. While quasars and BL Lacs 
generally exhibit a core fraction (defined Rc — Sc/ S total) Rc ~80%, the 30 galaxies in our 
sample show no tendency toward either high or low core fraction. The tendency towards 
higher core fraction in quasar and BL Lacs is expected from unified schemes since the core 
fraction is a strong function of orientation. The Kolmogorov-Smirnov (K-S) test disproves 
the null hypothesis that the distribution of quasars {dq) and the distribution of BL Lacs 
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(^bl) are the same as the distribution of galaxy core fractions {da)- When comparing dq 
to do the K-S test outputs the maximum value of the absolute difference between the two 
cumulative distribution functions (c.d.f.'s), D=0A5 with a probability, p, that quasars and 
galaxies arc drawn from the same population of 7.7 x 10^^. Similarly, when comparing c/bl 
to dc, we find a probability of 2.9 x 10"'^. Comparing c/bl to rfq we find p=0.41, indicating 
that these samples could come from the same parent distribution. Since the cores and jets 
generally have different spectral indices, the core fraction will depend weekly on redshift. As 
an example, a source with Rc — 0.75 at a redshift of zero would shift to Rc — 0.83 at z—1 
and Rc — 0.87 at z—2, assuming a flat spectrum core and steep {a — —0.7) jets. In Fig. 2 
we have not attempted any redshift correction so this could account for the minor difference 
between the shape of the distribution for the quasars and BL Lac objects. 

The distribution of Rc in Fig. 2 for the galaxies appears somewhat bimodal. Of the 15 
sources with Rc < 0.5 more than half (9/15) are CSOs. No obvious properties are identified 
with the high Rc population, although our knowledge of these sources is not complete. 

One might be concerned that this comparison uses a sample of quasars four to five times 
as large as the samples of BL Lacs and radio galaxies, respectively. This should not be a 
factor if we consider the effective number of data points A^e as defined by Press et al. (1992). 
Ne is defined as 



K^^^ (1) 

where Ni are the number of data points in the respective sets. As long as Ne is greater 
than or equal to four the probability estimate should be accurate. For the quasar— galaxy 
comparison Ne — 23.4, quasar— BL Lac Ne — 16.8, and for BL Lac— galaxy Ne — 12.0. 



4.2. Core and Jet Fractional Polarizations 



Fig. 3 shows the distribution of average fractional polarization in the cores of quasars, 
galaxies and BL Lac objects. From this distribution we see that quasars have typical core 
polarizations of a few percent, and most (87%) are detected. Galaxies are clearly less polar- 
ized than quasars on the parsec scale with only 3 of 30 (10%) of the sources found to have 
detectable polarization in their cores. No detectable polarized flux was found from any of 
the 9 CSOs observed. This is consistent with results from a VLBl follow-up survey of CSOs 
by Peck & Taylor (2000) who found less than 1.2 mJy of polarized flux from each of 21 
sources observed at 8.4 GHz with the VLBA. Generally we find less than 5 mJy of polarized 
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flux density and fractional polarizations of <0.5%. For a few CSOs with weak cores, our 
upper limits on the fractional polarization are quite high. In contrast the BL Lacs appear 
to have slightly more strongly polarized cores than quasars (K-S test gives p—0.20 that they 
are drawn from the same population). 

The K-S test only compares detections, and ignores the upper limits. Survival analysis 
(Feigelson & Nelson 1985) takes into account the upper limits in the comparison of two 
samples. Using the Penn State ASURV 1.2 package (Lavalley, Isobe, & Feigelson 1992) 
we tested the hypothesis that the BL Lac and quasar core fractional polarization results 
were drawn from the same parent distribution. Including the 15 upper limit quasar results 
decreased the likelihood that the two samples were drawn from the same parent distribution. 
The five two sample tests in ASURV 1.2 give an average probability of 1.5% that the null 
hypothesis is correct. This is significantly different from the K-S test result, and we find 
that this result is not affected by eliminating the one rather high upper limit quasar core 
fractional polarization of 2.7%. 

Fig. 4 shows the same set of distributions, but for the average fractional polarizations 
of the jet components. In this case the fractional polarization of the BL Lacs is significantly 
greater from that of the quasars (p=0.0037 that they are drawn from the same population). 
When the jet polarization is detected, sources have higher fractional polarization in the jets 
than in the cores. Even the galaxies have an average fractional polarization in the jets of 
10%, compared to typical fractional core polarizations of <1%. In Fig. 5 we show the images 
for all six radio galaxies with some detected polarization. The measurements range from 
marginal in 0600+442 and 0847+379 to extremely detailed in the bright, one-sided jet of 
1807+698. 

4.3. Compciring Alignments Between Jet Axis and EVPA 

Many claims have been made in the literature regarding comparisons of the orientation 
of the jet axis (6) with the EVPA of the jet (Xjct) oi core (Xcorc)- Below we discuss various 
tests to look for preferred alignments. Ideally these tests should be made after having 
removed any Faraday rotation due to the ISM of our Galaxy, or the environment local to 
the sources. Unfortunately, for these single-frequency observations it was not possible to 
estimate the Rotation Measure (RM). If significant RMs exist towards any of the source 
components (e.g., as found to be common in quasar cores by Taylor 1998, 2000) then that 
will tend to smear out any intrinsic correlation. 

The rotation measure of the cores of BL Lacs are typically 2—300 rad m~^, and those 
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of quasars range from 200 rad to several 1000 rad (Zavala & Taylor 2003). At 5 
GHz a rotation measure of 277 rad m~^ is enough to cause a turn of one radian. This would 
be sufficient to smear out any correlation in EVPA. The RMs in the jets of 16 BL Lacs, 
radio galaxies and quasars, with one exception, are all approximately 200 rad m~^, and thus 
would also smear out any correlation in ^ — xjet (Zavala & Taylor 2003). 

In Fig. 6 we present the distribution of ^ — Xcorc ^or quasars, galaxies and BL Lac objects. 
For the quasars we see that the distribution is not quite flat, and hints at an excess of sources 
near 80 degrees. These nearly misaligned EVPAs imply a magnetic flcid nearly perpendicular 
to the jet axis assuming optically thick synchrotron emission. The K-S test indicates only a 
0.3% chance that the distribution is the same as a flat distribution, suggesting the possibility 
of an intrinsic correlation between EVPA and 9. There is a somewhat greater probability 
(3.7%) that the quasars and BL Lacs are drawn from the same population. However, the 
BL Lac distribution itself is not significantly different from a flat distribution (17% chance 
that it is the same as flat) . 

To further investigate the nature of the misafigned core EVPAs we divided the quasars 
into those with short (length < 6 mas) jets and long (length > 6 mas) jets. From Fig. 7 we see 
that a jet length of 6 mas divides the quasars into two groups of roughly equal sizes since there 
are jet length peaks at 4 and 9 mas. We used ASURV's correlation tests (Isobe, Feigelson, & 
Nelson 1986) to examine the correlation that appears in Fig. 8 for the quasars with short jets. 
Although it appears that short jets are correlated with misaligned core EVPAs, a Kendall's 
tau test gives a less than 1 sigma signiflcance of a correlation. Additionally, plots of \0~Xcore\ 
versus jet length show no obvious correlation, regardless of jet length. Although there is no 
evidence for a correlation of quasar core EVPA afignment with jet length the K— S test of 
quasars versus a uniform distribution indicates that quasars have a preference for misaligned 
core EVPAs. As one can see from Fig. 8 the cores associated with short jets appear more 
misaligned (probability of a flat distribution is 0.009%, probabifity that they have the same 
distribution as sources with long jets is 1.7%), with a clear excess near 90 degrees. We also 
divided the quasars into those with a high {Rc > 85%) or low {Rc < 85%) core fraction. 
Core fraction seems much less important as a discriminator of the misaligned core EVPA 
population. 

One can also look for any dependence on alignment with redshift. In Fig. 9 we plot 
the alignments of the jet and core EVPAs with jet axis against redshift. No trends with 
redshift are apparent for either the jets or the cores. Correlation tests performed using the 
ASURV package verifled that no correlation of either core or jet misalignment with redshift 
is present. This was true for fractional polarization in the cores and jets as well. 

In Fig. 10 we present the distribution of 6* — Xjet for quasars, galaxies and BL Lac 
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objects for the jet components both near (< 6 mas) and far (> 6 mas) from the core. 
Contrary to claims in the hterature (Cawthorne et al. 1993; Homan et al. 2002) that quasars 
have preferred ahgnments with the magnetic field oriented parallel to the jet axis, we see 
no such correlation. (The probability of being drawn from a flat distribution is 34%.) We 
have also tried examining separately the near (< 6 mas) and far (> 6 mas) jet components 
and reach the same conclusions. Cawthorne et al. did apply a rotation measure correction, 
which might reveal a correlation our non-RM corrected results hide. However, their RMs 
were derived from the VLA, which does not correctly sample the parsec-scale RM of these 
AGN. Additionally, the RM corrections were primarily (12/17) less than 50 rad m~^, which 
amounts to a correction of 10° or less at 5 GHz. As shown in Zavala & Taylor (2003) these 
corrections were too small to account for the known parscc-scale RMs. Homan et al. did 
not apply a RM correction, but their observations were at 15 and 22 GHz. Their results are 
thus less susceptible to any rotation measure smearing. 

The distribution for BL Lac objects (Fig. 10) appears to have a lack of sources with 
EVPA perpendicular to the jet. From the K-S test we find an 8.0% chance of this population 
being drawn from a uniform 6 — Xjet distribution. A bimodal distribution was found by 
Gabuzda, Pushkarev & Cawthorne (2000) with peaks near 0° and 90°. Gabuzda et al. 
performed a similar analysis on a 1 Jy flux limited sample of BL Lac objects observed 
with VLBI polarimetry at 5 GHz. This is shown in Fig. 12 of Gabuzda, Pushkarev & 
Cawthorne (2000). A K-S test of the data in their Fig. 12 yields a 1.1% probability that the 
data are drawn from a uniform distribution. With assistance from Gabuzda (2002, private 
communication) we combined our BL Lac 9 — Xjet data with the 35 measurements from 
Gabuzda, Pushkarev & Cawthorne (2000) and tested this data against a uniform distribution 
in a K-S test. The combination of these two samples results in a 0.7% probability that the 
parent distribution is uniform. However, with just the 17 data points from our sample it is 
difficult to draw a conclusion on the nature of the underlying distribution, and wc cannot 
hope to isolate what type of source gives rise to the slight excess of sources with EVPA at 
small angles to the jet axis. 

In a 43 GHz polarimetry survey of fiat-spectrum sources from the Pearson-Readhead 
survey (Pearson & Rcadhead 1988), Lister (2001) found that the most strongly polarized 
quasar cores display EVPAs that are aligned with the jet axis. In Fig. 11 we examine 
the relationship between alignment and fractional polarization. No preferred alignment for 
the more polarized sources, as what Lister (2001) found, is evident in our sample. Also the 
preferred core EVPA orientation in our sample is for EVPAs perpendicular to the jet (Fig. 6). 
These two claims can be reconciled with the same projected magnetic field orientation if the 
cores are optically thin at 43 GHz and optically thick at 5 GHz. 
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Using the ASURV bivariate tests (Isobe, Feigelson, & Nelson 1986) we examined the 
hypothesis that there is a correlation between the core fractional polarization and the core 
EVPA alignment. The Kendall's tau test gives a probability of 82.7% and the Cox regression 
a probability of 96.9% that a correlation does not exist. Although there appears to be a 
slight excess of sources with misaligned core EVPAs at fractional polarizations greater than 
2.5%, this is not born out by the bivariate tests. The Cox regression gives a probability of 
94.7% and the Kendall's tau test a 43.6% probabihty that no correlation is present. The 
disparity probably results from small number statistics as only 20 sources have mcore > 2.5%. 

4.4. Fractional Polarization and Jet Length 

In Fig. 12 the core and average jet fractional polarization is plotted against the pro- 
jected jet length. The fractional jet polarization seems fairly smoothly distributed, but the 
fractional core polarization is anti-correlated with jet length. Sources with short jets seem 
to have higher core polarizations and there is a distinct lack of strongly polarized cores in 
sources with long jets. The censoring present in the lower panel of Fig. 12 led us to use 
survival analysis to test for a correlation for core EVPA alignment with jet length. Both the 
Kendall's tau and Spearman's rho tests give a weak (1.6 a) support to the anti-correlation 
which seems to be present. A possible explanation for this trend is that sources with long 
jets are viewed at angles farther from the linc-of-sight. Assuming the standard obscuring 
torus model, the angle to the core then traverses a higher- density region with magnetic 
fields causing Faraday depolarization. This interpretation is also consistent with the finding 
that the cores of BL Lacs (thought to be viewed at very small angles to the line-of-sight 
(Antonucci 1993)) are somewhat more strongly polarized than quasar cores. 

5. Cosmology Using the Size— Redshift Relation 

The 9 — z relation for compact radio sources has been studied in the past to place limits 
on go- Kellermann (1993) used a sample of 82 core-jet sources selected from the literature to 
find that the deceleration parameter is roughly 0.5, without appealing to source evolution. 
Gurvits, Kellermann, & Frey (1999) did a similar analysis using 330 5 GHz compact sources 
taken from the literature, finding that Qo — 0.21 ±0.30 if no linear size— luminosity, —redshift 
or —spectral index dependences are assumed. Gurvits, Kellermann, & Frey (1999) applied 
a regression model to 145 of the 330 sources with additional restrictions on luminosity and 
spectral index. With this subset of sources a deceleration parameter of go = 0.33 ± 0.11 was 
found. All of these results are consistent with Friedmann cosmologies with < go < 0.5 and 



-13- 



A=0. 

Figure 13 shows the O—z relationship for 103 core-jet sources from our complete sample 
and mirrors Fig. 1 in Kellermann (1993) and Fig. 5 in Gurvits, Kellermann, & Frey (1999). 
We measured angular sizes similarly to how Kellermann and Gurvits did, accounting for the 
decrease in brightness with redshift as described in § 3.3. CSOs have been excluded due to 

the proposed relation between their age and angular size, and BL Laccrtac objects have been 
excluded because of their proposed preferred inclination angles and to be consistent with the 
methods of Kellermann (1993). The redshifts range from 0.0172 < z < 3.469. With the 
exception of the slope between the first two redshift bins, our plot shows no clear dependence 
between angular size and redshift. Binning the data in redshift space rather than restricting 
each bin to equal numbers of sources does not greatly alter the plot. We found that plotting 
median values rather than mean values also does not change the plot's shape considerably. 
Thus contrary to previous claims our complete, flux-density limited sample, although not 
inconsistent with Friedmann cosmologies, cannot definitively rule out a steady state universe 
or place any limits on go- The few sources with higher angular sizes and reasonable error 
bars clustered around z=l.l are, however, suggestive of a general increase of angular size 
with redshift after z=l, which supports inflationary cosmologies. 

Although these compact radio sources with smaller redshifts and younger ages are not as 
subject to evolutionary effects, and although the physical parameters of their central engines 
arc thought to outweigh any effects due to possible variations of intcrgalactic medium with 
redshift, we find that these objects nevertheless cannot be used to place restrictions on the 
deceleration parameter and the geometry of our universe. The large error bars in Fig. 13 
prove these sources to be poor standard rods, possibly because of their rapid evolution 
compared to the Hubble time. It is possible that strong Doppler favoritism is skewing our 
results so that more objects with small inclinations to the line-of-sight and therefore greater, 
Doppler-boosted luminosities and smaller angular sizes are getting into the sample without 
being excluded by the fiux density requirement. Our neglect of Doppler boosting may also 
affect our measurement of angular size since we have defined object size by the contour line at 
2% of the peak in the total intensity map. Angular size might also be better measured taking 
into account jet bend, rather than assuming a one-dimensional jet in all cases. Lastly, it is 
possible that with a multi-frequency study using an angular size— spectral index relation, as 
is shown in Fig. 7 of Gurvits, Kellermann, & Frey (1999), the 9 — z correlation may become 
tighter. 
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5.1. Bent- Jet Sources and Other Interesting Properties 

Although wc did not do an exhaustive statistical analysis of the correlation between 
jet bending and other source properties, we do note two basic categories into which most 
bent sources fall. One class of bent jets typically shows little or no polarized intensity in the 
jet, while there may be polarized intensity at the source core. Some examples of this are 
0600+442, 0700+470, 0707+476, 0831+557, and 0843+575. The second and larger class of 
bent jets shows a tendency for the electric field to follow the curve of the bending axis so that 
the electric vector position angles are usually, but not always, perpendicular to the bending 
axis of the jet. Figure 14 shows 0133+479, 0627+532, 1151+408, 1459+480, 1619+491, and 
2351+456 which are all typical sources in this second category. A nice example of an E-field 
bending with and parallel to the jet axis is the quasar 0627+532. 

6. Conclusions and Future Work 

A new, fairly strong correlation between core EVPA and jet axis angle has been dis- 
covered in quasars with projected jet lengths shorter than 6 mas. These objects exhibit 
E-fields nearly perpendicular to the axis while those quasars with jet lengths longer than 
6 mas show no obvious 9 — Xcorc correlation. That this correlation exists in spite of the 
tendancy of quasar cores to have substantial Faraday RMs (Taylor 1998, 2000), implies that 
the actual correlation is even stronger than we observe. These findings suggest that quasars 
with axes anti-aligned with the E-field are not intrinsically different from those without this 
correlation; rather, the long-jet sources may be at larger angles to the hne-of-sight so that 
the angle to the core traverses a higher-density region with magnetic fields that produce 
sufficient Faraday rotation to smear out any intrinsic correlation. This reasoning also can 
explain the sparsely populated region of high projected jet lengths and high fractional core 
polarizations in Fig. 12. These results can be made consistent with the finding of Lister 
(2001), of electric vectors aligned with the jet axis at 43 GHz if the cores are optically thick 
at 5 GHz, but optically thin at 43 GHz. A magnetic field perpendicular to the jet direction 
could be produced by a strong transverse shock (e.g. Laing 1980). 

Unlike the preferred alignment found in quasar cores and contrary to previous claims in 
the literature (Cawthorne et al. 1993; Homan et al. 2002), quasar jets are not found to exhibit 
a, 9 — X relationship. This seems to be the case regardless of jet length, core dominance, or 
fractional polarization. 

In contrast to quasar cores, BL Lac cores in our sample show no strongly preferred 
alignment. Since BL Lacs are observed to have lower Faraday Rotation Measures than 
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quasars (Zavala & Taylor 2003), it is unlikely that the absence of correlation is due to 
Faraday rotation. BL Lacs may be intrinsically different from quasars (i.e., have a lower 
power), or their intrinsic high variability may make it impossible to find a correlation over 
just one epoch, with only a modest number of sources (17). Future multi-epoch analysis, 
and/or the analysis of a larger sample, is necessary to conclusively determine the existence of 
any preferred alignment. An interesting test to make in the future would be to use spectral 
information to look for preferential alignments in optically thick or thin cores. Contrary to 
claims made by Gabuzda, Pushkarev & Cawthorne (2000) of a bimodal 9 — Xjet distribution, 
the EVPA in BL Lac jets similarly show no correlation with the jet axis angle. 

We found no evolution of source properties (size, fractional polarization, or polarization 
angle) with redshift. This lack of evolution is reasonable if the physical parameters of the 
central engine outweigh any effects due to possible variations of intergalactic medium with 
redshift. 
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Fig. 1. — A selection of images from the sample showing a range of core- jet morphologies 
and the CSO 2352+495. Contours are drawn at — 4(T7, 4(T/, Sct/, 16ai, where ai is RMS of 
total intensity given in Table 1. Polarization vectors have lengths proportional to fractional 
polarization. Blanking of the polarization has been performed on pixels less than 4.6(7p in 
polarization, 9a"/ in total intensity, or 0.2% in fractional polarization. 
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Fig. 2. — Distribution of core fractions, Rc, in quasars, galaxies, BL Laccrtae objects and 
others, given as a percentage. Rc is the Stokes I flux density in the core divided by the 
integrated Stokes I flux density in the source, where we have calculated the core flux as 
described in § 3.3. The width of each bin is 4%. All sources are shown. 
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Fig. 3. — Distribution of the fractional polarization (vri) in the cores of quasars, galaxies, BL 
Lacertae objects, and others, given as a percentage. Upper limits are shown as dashed lines 
with the < symbol. The sohd line represents the sum of the upper hmits and detections 
in each bin. Upper limits were assumed when the polarized intensity ]) < 0.2% of the 
peak in the total intensity image, or when p < 3a, where a is the RMS noise calculated as 
described in § 3.3. The quasar distribution represents 106 sources with 15 upper limits. The 
galaxy distribution represents 30 sources with 27 upper limits. The BL Lacertae distribution 
represents 20 sources with no upper limits, and the distribution of other sources represents 
21 sources with 10 upper limits. The width of each bin is 1%. 
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Fig. 4. — Distribution of the average fractional polarization (m) in the jets of QuasarSj galax- 
ies, BL Lacertae objects, and others, given as a percentage. Only detected jet components 
have been averaged. (See § 3.3 for a definition of detected components.) The distribution 
of quasars represents 43 of the 106 we have in our sample. 21 quasars were excluded due 
to their naked core morphology, and 42 sources had no jet components with detected po- 
larizations. The distribution of galaxies represents 6 sources. Four were excluded as they 
were naked cores, and 20 had no detected jet components. The distribution of BL Lacertae 
objects represents 11 sources. Three BL Lacs were naked cores and 6 had no detected jet 
components. The distribution of others represents 7 sources. 6 were naked cores and 8 had 
undetected jet components. The width of each bin is 1%. 



-23- 



44 13 58.5i^ 0600+442 



O 



58.56 



58.54 



58.53 




06 04 35.631 35.630 35.629 35.628 
RIGHT ASCENSION (J2000) 



37 47 09.5(1^ 



09.48 



a 09.47 

a 



09.46 



0847+379 



09.45 \Q\ I I \ L 

08 50 24.732 24.731 24.730 24.729 

RIGHT ASCENSION (J2000) 



34 45 50.8(h 



3 50.79- 



O 
< 



u 
iij 
□ 



50.77 



50.76 j^ 
1013 49.616 



1010+350 




49.615 49.614 49.613 
RIGHT ASCENSION (J2000) 



49.612 



39 14 22.7(1- 



S 22.69 



O 
lu 
Q 



22.67 



22.66 



1223+395 




12 25 50.572 50.571 50.570 50.569 50.£ 

RIGHT ASCENSION (J2000) 



37 20 21 .6i^ 1456+375 



U 
lU 
Q 




14 58 44.796 44.795 44.794 44.793 

RIGHT ASCENSION (J2000) 



}49 28.1'h 



o 

28.10 



1807+698 




28.09 



18 06 50.684 50.682 50.680 50.678 50.676 50.674 50.672 
RIGHT ASCENSION (J2000) 



Fig. 5. — A plot of all six radio galaxies within our sample having some detected polarization. 
Contours are drawn at —4(7/, 4(7/, 8a i, 16a i, where ar is RMS of total intensity given in 
Table 1. Polarization vectors have lengths proportional to fractional polarization. Blanking 
of the polarization has been performed on pixels less than Sap in polarization, Aaj in total 
intensity, or 0.2% in fractional polarization. 
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Fig. 6. — Distribution of the difference between the jet axis angle, 9, and the electric vector 
position angle at the core, Xcore given in degrees. The distribution of quasars represents 
72 sources. 21 quasars were excluded due to their naked core morphologies, and 13 were 
excluded as they did not meet the detection criteria described in § 3.3. The distribution of 
galaxies represents just 3 sources, where 4 galaxies were found to have naked core morpholo- 
gies and 23 had undetected cores. 17 BL Lacertae objects are shown. Three naked core BL 
Lacs were excluded. The distribution of others represents 6 sources. 9 had undetected core 
polarizations, and 6 were naked cores. The width of each bin is 3°. 
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Fig. 7. — Distribution of uncorrected jet lengths for the 72 quasars shown in Fig. 6. The 
width of each bin is 1 mas. 
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Fig. 8. — Distribution of the difference between the jet axis angle, ^, and the electric vector 
position angle at the core, Xcore, for quasars. Top left: Q — Xcore relation for those quasars 
with uncorrected lengths of less than 6 mas. 26 quasars are shown. Top right: relation 
for quasars with lengths greater than 6 mas. 46 quasars are shown. Bottom left: Q — Xcore 
relation for those quasars with core fractions greater than 0.85. 16 quasars are shown. 
Bottom right: relation for quasars with core fractions less than 0.85. 56 quasars are shown. 
Here we measure length as the angular distance from the core to the farthest jet component, 
irrespective of jet bend. The width of each bin is 4°. 
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Fig. 9. — Top: Plot of the average difference between the jet axis angle, 9, and the electric 
vector position angle in the detected jet components of quasars, Xjet, as a function of redshift. 
We have averaged the jet components within 6 mas of the core separately from those farther 
than 6 mas from the core, so that any quasar with detected jet components both farther 
and nearer than 6 mas from the core is represented twice. 24 quasars with detected jet 
components within 6 mas of the core are shown, and 17 quasars with detected jet components 
farther than 6 mas of the core are shown. Bottom: Plot of |^ — Xcorcl for those quasars with 
detected core polarizations as a function of redshift. 68 quasars are shown. 
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Fig. 11. — Top: Plot of the average difference between the jet axis angle, 9, and the electric 
vector position angle of all detected jet components nearer (o) and farther (x) than 6 mas 
from the core as a function of the average fractional polarization of the quasar jet. Only 
those quasars with detected jet polarizations are shown. 36 quasars are represented with the 
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Fig. 12. — Plots of the average fractional polarizations of jet components (top) and the 
fractional core polarization (bottom) vs. projected jet length. The projected jet lengths and 
upper limits for the fractional core polarizations are calculated as described in § 3.3. The 
top plot shows 43 quasars, 6 galaxies, 11 BL Lacs and 7 others. The bottom plot shows data 
for 106 quasars, 30 galaxies, 20 BL Lacs and 21 others. 
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Fig. 13. — Plot of the redshift vs. uncorrected mean angular size in mas, where the mean 
angular size is calculated as described in § 3.3. Our definition of mean angular size mirrors 
the definition given in Kellermann (1993). Only those sources with core-jet morphologies and 
known redshifts are included. Additionally, BL Lacertae objects have been omitted to give 
a total of 103 sources. Each point represents roughly equal numbers of sources. (Starting 
with the lowest redshift bin, the bins contain 10, 10, 10, 11, 10, 10, 11, 10, 10 and 11 sources, 
respectively.) We show la error bars, where (t^ is the standard deviation of the redshifts in 
each redshift bin, and ay is the standard deviation of the angular sizes in each redshift bin. 
The solid curves represent the theoretical dependence of angular size on redshift assuming 
a Friedmann cosmology with Ho=50 km s~^ Mpc~^ and a source with an angular size of 41 
parsecs. 
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Fig. 14. — A selection of images from the sample showing those sources with moderate to 
large bends and some detected polarization. Contours, polarization vectors, and blanking is 
the same as Fig. 1. 
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Measured Source Properties 
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2.0 


— 


— 


— 


0.13 


0.10 


1851+488 


A 


2 


Qn 


— 


248 


5.0 


144 


— 


— 


— 


— 


0.12 


0.13 


1908+484 


A 


2 


Oj 


60 


96 


<0.3 


— 


11.8 


— 


— 


— 


0.10 


0.10 


1910+375 


A 


2 


Qj 


176 


268 


11.7 


73 


8.1 


20.1 


— 


76 


0.13 


0.83 


1924+507 


A 


2 


Qj 





217 


1.7 


6 


9.3 


2.2 


36 


— 


0.11 


0.11 


1928+738 


C 


P 


Qj 


159 


1664 


15.9 


79 


20.2 


8.6 


80 


53 


0.66 


0.18 


1943+546 


A 


1 


Gc 


84 


14 


<0.6 


— 


25.5 


— 


— 


— 


0.24 


0.09 


1946+708 


A 


2 


Gc 


32 


29 


<0.3 


— 


18.1 


— 


— 


— 


0.31 


0.09 


1954+513 


C 


P 


Qj 


125 


580 


1.4 


45 


14.8 


— 


— 


— 


0.18 


0.12 


2021+614 


C 


P 


Qc 


41 


91 


<2.5 


— 


6.4 


— 


— 


— 


0.36 


0.18 


2054+611 


A 


2 


Oj 


165 


253 


3.0 


12 


7.1 


7.0 


62 


83 


0.12 


0.11 


2116+818 


A 


E 


Gj 


153 


72 


<0.3 




8.2 








0.12 


0.09 


2136+824 


A 


2 


Qj 


139 


117 


2.7 


128 


35.3 


5.1 


40 


40 


0.18 


0.12 


2138+389 


B 


2 


Qj 


90 


124 


<0.3 




6.0 








0.13 


0.11 


2200+420 


C 


P 


BLj 


162 


1098 


19.5 


170 


15.2 








0.35 


0.16 


2214+350 


A 


1 


Qj 


7 


346 


1.9 


173 


7.3 








0.12 


0.10 


2235+731 


A 


2 


Oj 


32 


287 


1.1 


106 


3.9 


8.0 


57 




0.11 


0.10 


2238+410 


A 


2 


Qj 


133 


320 


7.2 


11 


4.9 








0.11 


0.10 



-37- 



Source 


E 


P 


C 




I * 

^core 


Pcore 


Xcore 


L 


Pjet* 


N* 


p* 


0-/ 


(7p 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


2310+385 


A 


2 


Qj 


60 


202 


2.1 


102 


4.8 


2.7 


53 




0.10 


0.11 


2319+444 


A 


2 


Qn 




303 


15.7 


134 




2.7 






0.11 


0.10 


2346+385 


A 


2 


Qj 


149 


344 


10.2 


51 


5.2 








0.14 


0.10 


2351+456 


C 


P 


Qj 


108 


768 


14.8 


66 


9.8 


4.5 


14 




0.25 


0.18 


2352+495 


A 


P 


Gc 


2 


8 


<0.9 




25.5 








0.33 


0.09 


2356+385 


A 


2 


Qn 




329 


3.6 


24 










0.12 


0.10 


2356+390 


A 


2 




51 


125 


8.5 


161 


10.8 


7.5 





28 


0.09 


0.09 



*A more complete table with all measured source properties, including jet components, can be found at 
www . aoc . nrao . edu/~gtaylor/c j f tab . text. 

Note. — Col. (1): B1950 source name. Col. (2): E, epoch - A:1998; B:1999; C:2000. Col. (3): P, 
parent sample - P:Pearson & Readhead 1988; 1: first Caltech-Jodrell Bank survey (CJl, Polatidis et al. 
1995); 2: second Caltcch-Jodrcll Bank survey (CJ2, Taylor et al. 1994); E: 18 sources not included in PR, 
CJl or CJ2 that complete the Caltech-Jodrell Bank flat-spectrum survey (CJF). Col. (4): C, class. Key to 
identifications: Q:quasar; G:galaxy; BL:BL Lac object; 0:other. Key to subscripts: j:core-jet morphology; 
n:naked core morphology; c:Compact Symmetric Object. Col. (5): jet axis angle (deg). Col. (6): total 
intensity at the core (mJy). Col. (7): total polarized intensity at the core (mJy); cores with detected 
polarized intensities have IcTtypicai ~ 0.2 mJy. Col. (8): electric vector position angle at the core (deg); 
cores that meet detection criteria have Icrtypicai ~ 4°. Col. (9): L,jet length (mas). Col. (10): average 
polarized intensity of detected jet components (mJy). Col. (11): N, average difference between jet axis angle 
and electric vector position angle in detected jet components within 6 mas of the core (deg). Col. (12): 
F, average difference between jet axis angle and electric vector position angle in detected jet components 
farther than 6 mas from the core (deg). Col. (13): total intensity RMS noise (mJy). Col. (14): polarized 
intensity RMS noise (mJy). 



